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Coherent transport: Landauer-Buttiker method 

Quantum junctions • Landauer formula • Multi-channel scattering and transport 
Multi-terminal systems • Problems 



The Landauer-Buttiker (LB) method (also known as the scattering method) estab- 
lishes the relation between the wave functions (scattering amplitudes) of electrons in 
a nanojunction and the conducting properties of this junction. The LB method can 
be applied to find the current through a noninteracting junction or through a junction 
with effectively noninteracting quasiparticles, for example if the mean-field descrip- 
tion is valid and the inelastic scattering is not essential. Such type of an electron 
transport is called coherent, because there is no phase-breaking and quantum inter- 
ference is preserved during the electron motion across the system. Actually, the LB 
method is now routinely applied to basic transport calculations through nanostruc- 
tures and single molecules. Besides, it is directly applicable in many semiconductor 
quantum dot systems with weak electron-electron interaction. Due to simplicity and 
generality of the LB method, it is now widely accepted and is in the base of our 
understanding of the coherent transport. 

In this chapter we consider the foundations and some applications of the LB 
method, related directly to the problems of quantum transport at nanoscale. At the 
beginning (section 1.1) we discuss a basic quantum junction model and introduce 
necessary results from the quantum scattering theory. Then we obtain the Landauer 
formula in the single-channel (and independent-channel) case and discuss its phys- 
ical sense (section 1.2). In section 1.3 the general multi-channel formalism is dis- 
cussed, and in section 1 .4 the multi-terminal case is considered. 

Some additional topics, such as noise, thermal transport, hybrid systems with 
superconducting and magnetic electrodes, are briefly discussed in section ??. Many 
other extensions of the method are beyond our consideration. One can mention the 
LB method in combination with random matrix theory, localization theory, the at- 
tempts to introduce the analogous scattering description for interacting and dissipa- 
tive systems. However, in spite of such popularity of the LB method, one should 
remember, that in its canonical form it is applicable only for noninteracting systems. 
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1.1 Quantum junctions 
1.1.1 Electrodes, leads, scatterer 

As we already discussed in the Introduction, a nanojunction consists of three parts: 
left and right large equilibrium electrodes and central nanoscale region. In this chap- 
ter we neglect possible discrete nature of the central region and consider the nano- 
junction as a continuous medium. Moreover, we assume, that transport is completely 
coherent between the electrodes. It means that we neglect all inelastic effects in- 
side the central region. The electrodes, however, are equilibrium and incoherent by 
definition, because the equilibrium state is established as a result of energy relax- 
ation due to inelastic scattering. If an electron comes from the central region into 
the electrode, it is thermalized and any phase information is lost. On the other hand, 
electrons, which enter to the central region from the electrodes with random phases, 
keep the phase information until they return back to the electrodes, but their distri- 
bution function in the central part is determined by the boundary conditions and is 
nonequilibrium at finite voltage. For further consideration it is convenient to divide 
the central region into quantum leads and scattering region, or scatterer (Fig. 1.1). 

Transport through a coherent region can be described by a wave function. The 
leads serve as quantum waveguides for electrons. They connect the electrodes with 
the scattering region and are assumed to have a well defined mode structure: incom- 
ing (from the left or right, s = L,R, electrode to the scatterer) modes y/ s +(r) and 
outgoing modes y/ s - (r) can be defined. In the simplest case the leads are noninter- 
acting, but nonequilibrium at finite biased voltage. The scattering region can be as 
simple as one tunneling barrier or impurity, or as complex as an interacting nonequi- 
librium molecule. In the continuous case the scatterer is described by some potential 
U (r), it results from impurity scattering or external potential. 

Note, that the boundaries between electrodes, leads and scatterer are to a certain 
extent arbitrary. It is important only to take all conduction modes (transport channels) 
into account. 



Fig. 1.1. Schematic picture of a quantum junction: the scattering region is connected to the 
electrodes trough the quantum leads. 
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1.1.2 Transport channels 

Let us introduce the important concept of transport channels. To start, we simplify 
the problem, assuming that the motion of electrons is effectively one-dimensional. 
For this purpose we factorize the wave function with some energy E in the following 
form: 

y £ (r)=2>„(x,y,z)y/„ £ (z), (1.1) 

n 

where z axis denote the direction of electron motion, <j> n (x,y,z) describes the trans- 
verse structure, n is called mode or channel index. At the first glance, we only com- 
plicated the problem, because instead of one unknown function we got many. But 
the trick is that the solution for (j) n (x,y,z) can be easily obtained in many cases and 
we should solve only the equation for ^r nE {z). For example, in 3D layered systems a 
corresponding solution is 

«Mr)=£e^^Wz), (1-2) 

k 

where k is the wave vector parallel to the layers. Initial 3D problem is reduced in 
this case to a ID problem for the function y/^z). Similar situation takes place in 
the effectively ID or 2D spatially quantized systems (electronic waveguides), where 
<i>nE{x,y,z) describes different transverse modes at different n. 

To understand better, how transport channels are formed, let us consider two ex- 
amples (Fig. 1.2): a constant cross-section waveguide and an adiabatic waveguide 
(quantum point contact). Such quantum waveguides for electrons can be three- 
dimensional or formed from 2D electron gas, in the last case the motion in y direction 
is quantized and y coordinate is absent in the equations below. 

The wave functions are calculated from the Schrodinger equation 



h 2 ( d 2 . d 2 d 2 \ 



U(x,y,z) \ W E {x,y,z)=EW E {x,y,z), (1.3) 



2m \dx 2 dy 2 dz 2 ) 
where U(x,y,z) is the external potential, defining the geometry of the junction. 
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For waveguides with constant cross-section U(x,y) is z-independent, and the so- 
lution can be presented in the form (1.1) exactly as 

f £ (r)=£0„(x, y y^, (1.4) 

n 

with the wave vectors hk„ — ±y/2m(E — E„), E n are the transverse energy eigenval- 
ues. From this expression it is clear, that the number of propagating modes (called 
open channels) at given energy E equals to the number of transverse modes with 
E < E: 

N ch (E)=Y,0(E-E n ). (1.5) 

n 

AtE„>E the wave number is imaginary and this mode can exist only at the end of 
the waveguide, these channels are called closed. It is important that N c h is a function 
of energy and can be changed if the Fermi energy or the gate voltage are changed. 



Adiabatic junction 

Now we consider a waveguide with the variable cross-section. In the adiabatic ap- 
proximation the spatial variation in the z-direction is much slower than in transverse 
directions. Thus, the second derivative with respect to z is small and it is reasonable 
to consider the Scrodinger equation for the transverse eigenfunctions <j> n {x,y) with z 
as a parameter 



2^m {lL? + + U ( x ^^ z )^n{x,y 1 z) =E n (z)<j) n (x,y,z), (1.6) 

E n (z) is the (position-dependent) energy eigenvalue. If we substitute the Ansatz (1.1) 
with these eigenfunctions into (1.3) and neglect again the derivatives of <j) n (x,y,z) in 
the z-direction, we get the ID equation for the function y/„£(z): 

~^+E n {z)^\lf nE {z)=E\if nE (z), (1.7) 

where the energy E n (z) plays the role of an effective potential. 

Consider a simple example, namely a 2D junction with the hard-wall potential. 
Assume, that the size of the waveguide in x-direction is d(z). The transverse solution 
(at fixed z) is 

^ z)= \lW) sm \ — W) — J' (L8) 

and corresponding energy eigenvalue, e.g. the effective potential for z motion is 

h 2 7l 2 n 2 

E ^ = ^k Y <L9) 

If we plot now this energy as a function of z (Fig- 1.3), we can determine, which 
channels are open at given energy E. It will be all channels with the maximum E n (z) 
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Fig. 1.3. Effective potential energy E n (z). The open transport channels are shown by solid 
lines, closed one - by dashed lines. The dotted line shows d(z). 



being smaller than E. Note, that if the maximum E n {z) is close to the energy, the 
channel is only partially open, because both scattering and transmission are possible 
with some probabilities. For adiabatic junctions, however, this energy interval is very 
small, and all channels can be considered to be open or closed. 

The condition of adiabaticity in two dimensional junctions can be written as 
dd{z)/dz <C Xf/cI{z). On the other hand, the position-dependent number of open 
channels is N(z) — k F d(z) jit and one gets 



dd(z) 



1 



dz N(z) ' 



(1.10) 



This condition shows, that it is much easier to obtain the adiabatic transport for nar- 
row junctions with small number of channels. 

In the above considered examples the channels are independent, there are no 
transitions between different modes. This type of transport can be named mode- 
conserving. Below in this and next sections we consider only single channel and 
independent channels. More general situation with inter-mode scattering will be con- 
sider in the section 1.3. 



1.1.3 Reflection and transmission 

As one can see, the main formal problem to be solved in the Landauer theory is 
the single-particle scattering problem. The conductance is determined then from the 
elements of the transmission matrix. For this reason, the LB method is called also 
a scattering method. In the following sections we introduce the basic notions of the 
quantum scattering theory, and formulate them in the form convenient for further 
calculations. 

Now consider a one-dimensional potential (schematic in Fig. 1.4) which is con- 
stant far from the scattering region: U(z^ — °°) = 0, U(z —¥ °°) = Uq. Thus, the wave 
function in the state with the energy E z (the full energy E =E n +E Z is assumed to be 
the sum of E z and the transverse energy E n ), and incident from the left, is 
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Fig. 1.4. One-dimensional potential. Reflection and transmission. 



y( z ^ _«) =A + e'' fc +A_e-' fa , (1.11) 
y( z ^+oo)=JJ_e a '* 1 (1.12) 

where A + is the incoming and A_, Z? are the outgoing waves (Fig. 1.4). Wave vec- 
tors k and k! are defined as 

s/2mE z 

k= — - — , (1.13) 
n 

k , = v — v* (v (114) 

n 

To determine the transmission and reflection probabilities, one should use the 
conservation of probability flux density ("current") 

ih 

f = ^ [v(r)Vy/» - y/*(r)Vy/(r)] . (1.15) 
From the Schrodinger equation the conservation law follows: 

«^+divf = 0, (1.16) 
at 

which in the stationary case is simply 

divf=0. (1.17) 

For the plane wave Ae' kz the probability flux density ("current") is reduced to 

hk 

f=-\A\ 2 = v\A\ 2 , (1.18) 
m 

where we introduced the velocity v = hk/m. 

Now we are ready to define the transmission coefficient as the ratio of the trans- 
mitted to incident probability flux 



ftran = V°q\B-\ 2 = V\R- 
finc V^\A + \2 k |A_| 



T(E Z ) = J Jpi= y"-\ 2 = ^fPz, (1.19) 
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Fig. 1.5. Single barrier. General scattering problem. 



and the reflection coefficient as the ratio of the reflected to the incident probability 
flux 

fref _ \A-\ 2 



R(E Z ) 



fine |A H 



From the probability flux conservation it follows that 

T(E Z )+R(E Z ) = 1. 



(1.20) 



(1.21) 



Note, that if we consider the incident wave from the right side of the barrier, the 
transmission coefficient is the same at the same energy 



(1.22) 



1.1.4 Scattering and transfer matrices 

Now we consider the general scattering (or transmission) problem, assuming that 
there are incoming modes from the left and from the right sides of the barrier 
(Fig. 1.5). 



8 -potential 

We start from the 5-potential (some other barriers are considered in Problems 1.2-4) 

U(z) = aS(z). (1.23) 

The solution is given by 



, , A+e ikz +A-e-' kz , z < 0, 
VW-\ B _ e ikz +B e -fc )Z>0) 



(1.24) 



J2mE- 

where A+, B + are incoming and A_, B are outgoing waves, k = . The 

« 

boundary conditions at the 5-potential are 

y/(0-) = y/(0+), 



y/(0+)-y/(0-) = ^y/(0). 



(1.25) 
(1.26) 
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. / h 2 k\ 

We can represent this boundary condition using the transfer matrix M: I K = J 



M 



K K 
- 1 - 



5_ 



(1.27) 



or, alternatively, the scattering matrix S: 
= 51 



1 iK 

iK-1 iK-1 \ ( A A 

iK ___]_ I \B + 

iK-1 iK-1 



(1.28) 



The transfer matrix relates the states from two sides of the barrier, while the scatter- 
ing matrix relates the amplitudes of outgoing waves to the amplitudes of incoming 
waves. 

To find the transmission and reflection coefficients we set now B + = 0, then 



T(E) = 





IB- 1 2 




1 


'v L 




|A + | 2 




|M n |2 - 


vr. 



R(E) 



\A + \ 2 



\M 2 i\ 
" \M n \ 
T+R = l. 



\Su 



1^21 1 2 = 



1 



K 2 



1+K 2 ' 



1+K 2 ' 



(1.29) 
(1.30) 



We use here the square brackets to show the terms, which appear for the junctions 
with different left and right velocities v, see the Problems 1.3, 1.4 for detailed calcu- 
lations. 

The descriptions in terms of scattering or transfer matrices are completely equiv- 
alent and the choice is only dependent on the convenience and the problem to be 
solved. Typically, in the end of a calculation the 5-matrix should be obtained, be- 
cause it determines the conductivity by the Landauer formula. But to calculate the 
scattering by the complex system, the M-matrix can be convenient, as we shall see 
in the next section. 



Some general properties of S and M matrices 

In general case (asymmetric junction, magnetic field, etc.) the scattering matrix is 
usually represented as 

'rt' 
t r> 



S = 



(1-31) 



The reflection and transmission coefficients from left to right are 

T L ^ R (E) = ^\t\ 2 ; R L ^ R (E) = \r\ 2 , 
Vr 



(1.32) 



and from right to left 



W£) = JV| 2 ; R R ^ L (E) = \r'\ 2 
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The transfer matrix can be written explicitly in terms of transmission and reflec- 
tion amplitudes r, r', t, t' as 



M 



1 _/ 

r tr - rr 



(1.34) 



t t 

as well the scattering matrix in terms of the elements of M as 

/ M 2 \ M n M 2 2-M n M 2 \\ 

S = 



1 



(1.35) 



For symmetric leads (Vl = Vr) the S-matrix is unitary 

SS f = I, 

so that 



kl 2 + kl 2 = i, 



\r'\ 2 + \t'\ 2 = l, r*t' + t*r' 



0. 



(1.36) 



(1.37) 



For asymmetric leads (Vl ^ Vr), however, the S-matrix (as we defined it up to 
now) is obviously not unitary, because of additional velocity factors in the expression 
(1.19) for the transmission coefficient. If one needs to have the unitary matrix, it is 
possible to redefine it as 



(1.38) 



Or, equivalently, redefine the amplitudes of plane waves as Ae lkz = (A' / y/v)e' kz , and 
the flux density 



f=-\A\ 2 = v\A\ 2 = \A 



'|2 



(1-39) 



Below we always assume that the S-matrix is unitary, and the reflection and trans- 
mission coefficients acquire the simple form 



T L ^ R (E) = \t\ 2 ; R L ^ R (E) = \r\ 



(E) 



\t'\ 2 ; R R ^ L (E) 



J\2 



(1.40) 



In the time-reversal case (no magnetic field) the product of the S-matrix with its 
complex conjugate is the unity matrix 



SS* = I, 



(1.41) 



from which follows 



30 1 Coherent transport: Landauer-Buttiker method 







b. a: 

— ► — 




B'_ 


A. 


''i, 


B A' 






■* — 


■ t'> 















z 



Fig. 1.6. Double barrier. 



\rf = \r\ 2 , |r| 2 +rY = l, \tf = \t\ 2 , (1.42) 

and 

R L ^ R {E) = R R ^ L (E); T L ^ R (E) = T R ^ L (E). (1.43) 

Note, finally, that an unitary matrix in the time reversal case is symmetric, and 
|r| = |r'|, so that the S-matrix can be represented in the form 




(1.44) 



A series of scatterers: transfer matrix method 

If one has a series of scatterers, the calculation of the resulting transfer matrix can 
be simplified by the use of the transfer matrices of single scatteres. Consider two 
sequential barriers with transfer matrices M and M (Fig. 1 .6), so that 

In the case when both M and M are calculated independently, assuming that the 
barrier is placed at z — 0, the outgoing coefficients B and incoming coefficients A' 
are related by the propagation transfer matrix Ml 

(t)=M£)=Co^)C£)- 

To show that in the most simple way, note that B-e lkz and A' + e tkz describe the 
same plane wave in two different points z — and z — L, the phase difference is 
obviously kL. 

Finally we can write 

(l:)=^(t) = ^'(t)- (L47) 

Thus, the transfer matrix for a sequence of barriers can be defined as the product 
of the particular transfer matrices of the barriers and the propagating transfer matrices 

M T =M x M\ 2 M 2 ..M n M^ n+l M n+l ..M N ^M^_ lN M N . (1.48) 
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Breit-Wigner formula 

As an example, we can calculate the transmission coefficient for the double-barrier 
system. We need only Mm because it determines T (see (1.29)) and R = 1 — T. The 
transfer matrix Mr for a two-barrier structure is 

M TU =M n M' n e~ ikL +Mi 2 M! 2l e ikL (1.50) 
In the case of two identical barriers, for transmission coefficient we find 

t ^ = t 2 + ^Ul-9Y (L51) 

where 9 is the phase of the complex M\\. T\ and R\ are transmission and reflection 
coefficients of the single barrier. 

From this general expression one can see the important property of two-barrier 
structures: there are transmission resonances, at some specific energies E n the trans- 
mission coefficient is large (T(E„) = 1 in symmetric structures), while between res- 
onances it can be small. 

i 1 ma 

When the barriers are o -functions M\\ = 1 H — , 9 = arctan — = arctan and 

K' K h 2 k 

the equation for resonances (T = 1) is 

h 2 k 

tan/fcL = . (1.52) 

ma 

Close to the resonance, around one of the resonance energies E nj the transmission 
coefficient has a Lorentzian form 

w- wm* ' ,L53) 



where the width r n is given for two 5-barriers as 

mL 2 R\ 



, 2h 2 E n T 2 \ l/2 
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1.2 Landauer formula 
1.2.1 Single-channel formulas 

The main idea of the scattering approach to the conductance was first formulated 
by Rolf Landauer [1, 2]. He proposed, that the conductance of some segment of a 
ID channel with elastic scatterers is determined by the quantum mechanical prob- 
abilities of transmission (T) and reflection (R = 1 — T) through this segment. It 
should be noted, that Landauer considered the local resistance of a system (the zero- 
temperature residual resistance), but not the resistance of a quantum system between 
two equilibrium electrodes. As a result, he got for the zero temperature one-channel 
(effectively one-dimensional) conductance the so-called "first Landauer formula" 

The result, which seems to be reasonable at least in two limiting cases. At small 
transmission T — > 0, the conductance is also small and proportional to T, the result, 
which is well known from the perturbation theory. In the opposite case, when T — >• 1, 
S->0, there is no scattering at all, so that the conductance should go to infinity, in 
agreement with (1.55). To take into account the spin degeneracy in this formula, one 
has to multiply the conductance (1.55) by 2. 

However, the further investigations [3, 4] show that the conductance of a ID 
system, calculated by the exact linear response method, can have also quite different 
form (depending on the boundary conditions) 

e 2 

G=—T. (1.56) 
h 

This conductance is finite even in the case of the perfectly transmitted junction 
(T = 1). Actually, there is no contradiction between these two formulas. It was shown 
that both are reasonable and give the same current, but correspond to the voltages, 
defined between different points. As we shall see below, the key difference between 
the formulas (1.55) and (1.56) is that the first one is for the conductance inside the 
junction (between points A and B, see Fig. 1.9 below), while the second gives the 
conductane related to the equilibrium electrodes (between points L and R in the 
Fig. 1.9). In the section 1.2.4 we obtain both formulas and discuss the relation be- 
tween them. The puzzle with finite resistance at T — > 1 is also understood, it is clear 
now that the current through a junction is always accompanied by the voltage drop 
at the boundaries between electrodes and leads. The physical reason is that the num- 
ber of open electron transport channels is limited, while many other electrons reflect 
from the junction and create some charge distribution. Not so obvious is, however, 
that this contact resistance has the universal value R c =h/e for one spinless channel. 

For the transport problems, considered in this book, the second type of the Lan- 
dauer formula is more important usually. Besides, the first type formulas are not exact 
for finite-size nanostructures, because they are dependent on the particular electrical 
potential distribution inside the junction. 
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The important question, discussed in connection with the Landauer resistance, is 
the origin of dissipation in this approach. Indeed, finite dc current at finite dc voltage 
means that the energy is permanently dissipated. On the other hand, we consider only 
elastic scattering, so that the energy can not be dissipated in the scattering process. 
This problem is closely related to the phenomena of the residual resistance at low 
temperature, caused by impurities. In both cases we should introduce some thermal- 
isation. In the case of transport between the equilibrium electrodes, this problem is 
resolved quite easy, the energy is dissipated in the electrodes, the details of the dis- 
sipation are not relevant. More precisely, the incoming from the electrodes particles 
are equilibrium distributed, while outgoing particles propagate into the electrodes 
and are thermalised here. 

At finite temperature and finite voltage the Landauer formula (1.56) is trans- 
formed into the more general formula for the current: 

e 1'°° 

I(V) = -J_J(E,V)[ME + eV)-f (E)]dE, (1.57) 
where f°(E) is the Fermi-Dirac distribution function with the Fermi energy £/r: 



ME) = 




(1.58) 



1.2.2 Heuristic derivation 

Now we are ready to see in detail, how the transmission coefficient can be used to 
calculate the current through a quantum junction. We start from the mode-conserving 
scattering and use here the heuristic arguments. More rigorous methods are summa- 
rized in section 1.3. 

From the scattering picture it follows that all particles, coming from the left elec- 
trode, are transmitted through the junction with the probability T(n,k z ) and, after 
that, their excess energy, phase coherence, and the memory of their previous state 
are lost in the right electrode. We assume in all cases, that an electron can go without 
scattering from the lead into the electrode, thus for incoming from the left electrons 
there are only two possibilities: to go into right electrode with the probability T or to 
return back to the left electrode with the probability R. The same property takes place 
for all particles coming from the right and transmitted to the left. Transport through 
the junction is coherent in this model, energy E and transverse quantum number n 
are conserved (the case of the multi-channel scattering, when n is not conserved, 
will be considered later). Irreversibility is introduced through the relaxation in the 
electrodes. The main assumption is that the right-moving particles in the left lead 
are populated with the equilibrium distribution function of the left electrode f[ q (E) 
and the left-moving particles in the right lead are populated with the equilibrium 
distribution function of the right electrode f^(E) (see Fig. 1.7). 

According to this model, the current of electrons, which enter from the left elec- 
trode is determined by the following exspression 
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f°° dk 
J L ^R = eT T L ^ R (n,k z )v L (n,k z )f L {n,k z )-^, (1.59) 
„Jo 2k 

where v>L(n,k z ) is the group velocity of the particle with momentum k z , //,(«, is 
the distribution function, the form of this function is considered below. The inte- 
gration is only for right-moving particles with k z > 0. Note, that it is not necessary 
to multiply this expression additionally by the factor like (1 — f R ((n,k z )) as in the 
tunneling "golden rule" theory, because this factor describes the number of empty 
states in the right equilibrium electrode and should be included when the transition 
between left and right states is considered. Instead, in our approach we consider 
scattering states in the leads, which formally can be extended in the electrodes. The 
transmission coefficient from the left to the right is simply the probability to find a 
particle in the right part of this state. 
Taking into account that 

dE z {k z ) dE(k z ) 

ufe) = ^r = W (L60) 

where E(k z ) = E n +E z (k z ) is the full energy, we obtain 

Il^r=jY[ T L ^ R {n,E)f L (E)dE, (1.61) 
and similar expression for the current of right-incoming electrons 

7 *^ = ?E/ T R ^ L {n,E)f R (E)dE. (1.62) 

Note that the integration in this expressions is done from the bottom of conduction 
band E nL ^ R y Taking into account the symmetry of transmission coefficients (1.22) 
we get the expression for the current 

/ = 



~ h L J_Jn(E) [f L \E)-f R \E)]dE. (1.63) 




Fig. 1.7. Left-moving and right-moving particles in a perfect wire (left) and in a wire with 
scatterer (energy diagram). 
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The limits of integration over E can be taken infinite, because the closed channels 
have T„(E) = and do not contribute to the current. 

Finally, the distribution functions in this expression should be discussed. There 
are different possibilities to create a nonequilibrium state of the junction. In equi- 
librium the electro-chemical potential p, = ji + e(p should be the same in both elec- 
trodes. Here jl is the (internal) chemical potential, which determines the filling of 
electron bands in the electrodes, and <p is the electrostatic potential. One can create 
a difference of only (internal) chemical potentials (Fig. 1.8, left) if populates one of 
the electrodes with extra particles. This case, however, is quite difficult to realize in 
nanostructures, because any change of the particle density causes the change in the 
electric field. Moreover, typically the external voltage is applied to the electrodes, 
while the (internal) chemical potentials of the electrodes far from the junction are 
not changed, jii = jj,R = (Fig. 1.8, right). More generally, one can say that the 
difference in the electro-chemical potentials between two points taken inside the 
equilibrium electrodes, is always produced by the external voltage (fii — p,R = eV). 
To determine the exact distribution of the charge density and electrostatic potential 
near and inside the junction, the self-consistent solution of the coupled Schrodinger 
and Poisson equations is necessary. In this case the expression (1.63) should be used 
with care when the voltage is not small. Indeed, the potential U(z) is now a function 
of the applied voltage, and consequently the transmission coefficient is a function of 
the voltage too. 

The distribution functions in the general case are 



f L {E) 



1 



exp( J? -V^ )+ 1 



f R (E) = 



1 



exp [ E -^ T e fK 



(1.64) 



But usually the simplified form can be used, with explicitly written external voltage 
1{ y) = iL r UE,V)[fo(E+eV)-f (E)}dE. (1.65) 



Where f°{E) is the Fermi-Dirac distribution function with the equilibrium chemical 
potentials fi L = fi R = Sf: 



Hi 





U(z) 








U(z) 












■'} 



































Fig. 1.8. Energy diagrams for chemical potential difference (left picture, the number of elec- 
trons in the band is changed, but not the potential) and voltage difference (right picture, the 
electron band is shifted up, the potential is modified). 
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ME) = 




(1.66) 



The distribution functions in the electrodes are the functions of energy E only, 
thus one can introduce the transmission function 

T(E)=Y,T n (E). (1.67) 

n 

and obtain finally 

el""- 

I(V) = - J_J( E 'V) [fo(E + eV) -f Q (E)]dE. (1.68) 

This formula can be wrong, however, if the external magnetic field is applied, be- 
cause the magnetic field violates the time-reversal symmetry and the relation (1.22) 
may be violated too. 

The conductance at zero temperature is 

2 

G= e -rY. T n(£F). (1-69) 

" n 



1.2.3 Conductance quantization 
Perfect wire 

Consider now the conductance of a perfect wire adiabatically coupled to two elec- 
trodes. "Perfect wire" means that there are several open reflectionless channels with 
transmission coefficient T n (E) = 1. Thus all right going electrons inside the junction 
are populated only by the left electrode and left-going electrons are populated only 
by the right elecrode (Fig. 1 .7, left). We can say that right moving electrons have the 
(pseudo-) electro-chemical potential of the left elecrode fit, while left moving elec- 
trons of the right electrode Of course, the state of electrons inside the wire is not 
equilibrium, and these "left" and "right" chemical potentials give the number and en- 
ergy of corresponding particles in the channel, but they are not usual thermodynamic 
potentials. 

Now we simply use the expression for the current (1.68). The distribution func- 
tions in the electrodes at zero-temperature are the step-functions 

f L (E,V) = 8Qi+eV-E), (1.70) 
f R (E) = 9(n-E), (1.71) 

and for the current we obtain 

7 ( y ) = ?E / T n (E,V)[d(n + eV-E)-8(n-E)]dE = 

e ^ rH-E n +eV „1 
= zL T n (E,V)dE=-NV, (1.72) 
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where we used T„(E,V) = 1, and N is the number of open channels between p.L = 
fi+eV and p.R = fi. For the conductance one has 

e 2 

G=—N. (1.73) 
h 

It is accepted to call the conductance of a single-channel perfect wire with spin 
as the conductance quantum 

G ° = T~ ~ 11 M ^ S = 7 ' 748 ' 10_5 Q l ~ 12^00 Q (L?4) 
The corresponding resistance is 

R = ^2=12.9kQ. (1.75) 

Where does the resistance of a perfect wire come from? The origin of this resis- 
tance is in the mismatch between the large number of modes in the electrodes and a 
few channels in the wire. So this is not the resistance of a perfect wire, but rather the 
contact resistance of the interface between electrodes and wire. 



Quantum point contact 

In quantum point contacts (QPC), which have usually the adiabatic form, the con- 
ductance at low temperatures is quantized in accordance with (1.73). In the spin- 
degenerate case it can be written as, 

2e 2 

G=—"£e(E F -E n ), (1.76) 

" n 

where Ep is the Fermi energy, and E n is the maximum of the transverse energy E n (z). 
The Fermi energy in 2D electron gas can be changed by the gate voltage V g , in this 
way the conductance quantization was observed experimentally in the form of steps 
at the function G(V g ). 

At finite temperature the conductance steps are smeared. Besides, the steps are 
not perfect, if the junction is not adiabatic. This can be seen from the exactly solvable 
model with the potential 

V(x,z) = ^m(O 2 x 2 + V - l -m(o}z 2 . (1.77) 
The transmission coefficients have a simple form [5]: 

= exp [-2% (E - V - (n + 1 /2)h(O x ) /(h(O z )] + 1 ' (1 ' ?8) 



At co z -C co x we return to the adiabatic approximation and well defined steps. 
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Classical point contact 

It is interesting to compare the quantum conductance (1.73) with the conductance 
of a classical point contact with large width do >• Xp, known as Sharvin conduc- 
tance [6]. Following [7], this conductance for 2D ballistic channel with the width do 
between two Fermi gases can be written as 

et>F dn 

I=—doj-eV, (1.79) 

In 2D electron gas dn/d/j. = m/nh 2 , and we obtain (with spin degeneracy) 

G s = 2 4 k -^. (1.80) 
h % 

From quantum mechanical point of view Icpdo/n is the number of transverse chan- 
nels N. 



1.2.4 Contact resistance 

Consider now the single-channel case with the imperfect transmission T ^ 1, repeat- 
ing the same calculation as in (1.72) we obtain 

2 

I= e T T(fL L -fL R ) = e -TV 7 (1.81) 
n n 

e 2 

G=—T. (1.82) 
n 

This is the conductance between the reservoirs, e.g. between some two points "L" 
and "R" inside the electrodes (see Fig. 1.9). Now consider two other points "A" and 
"B" inside the leads. The distribution functions and corresponding "electro-chemical 
potentials" (these potentials are not true potentials, but give the correct number and 
energy of electrons, as we discussed before) are different for left and right moving 




Fig. 1.9. The points of voltage measurement: L,R in the equilibrium electrodes; A,B inside the 
leads. 
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electrons. Now, however, these potentials are different also at different sides of the 
scatterer (Fig. 1.7, right). The potential p^ of the right moving electrons is equal to 
pi only in the left part of the wire, as well as p^ = p R in the right part. All other 
electro-chemical potentials are modified by the reflection from the barrier. Assume, 
that one can approximate the charge redistribution in the leads due to scattering by 
some quasiequilibrium distributions with corrsponding pseudo-potentials p. For ex- 
ample, only the part of right moving electrons is transmitted through the barrier and 
corresponding potential should be TpL, but additionally (1 — T)p R are reflected and 
move back. Finally, we obtain 

p^=p L , p^ = Tp L + (l-T)p R , (1.83) 
p^=p R , Pt = Tp R + (l-T)p L . (1.84) 

The difference of both "left moving" and "right moving" chemical potentials 
across the barrier is the same 

P-L ~ P-R = P-L ~ P-R = (1-T)(Pl-Pr)- (1-85) 

We can identify this potential difference with the potential drop between points A 
andB 

eV AB = (\-T)(p L -p R ). (1.86) 
Thus we can define the conductance (with the current (1.81)) 

r>- 1 - el T - elj aw 

which is exactly "the first Landauer formula" (1.55). The voltage Vab appears as a 
result of charge redistribution around the scatterer. Not surprising that for perfect 
wire with T = 1 and R~0 this conductance is infinite. 

The conductances (1.82) and (1.87) obey the following relation: 

This result can be understood in the following way. G 1 can be considered as the full 
resistance of the junction, consisted from two sequential resistances of the scatterer 
(G'~ ! ) and of the contact resistance of the perfect wire (h/e 2 ). 

Consider additionally the conductance of the incoherent series of N scatterers, 
each having the transmission coefficient T\. If the phase coherence is broken, one 
should summarize the probabilities of transmission instead of the quantum ampli- 
tudes. Thus, the transfer matrix method does not work i this case. Instead we use 
the probability theory. Let us consider first only two scatterers with the transmission 
coefficients T\ and T 2 . The probability of transmission trough both scatterers T is 
calculated as the sum of all possible (re)scattering processes 

T = 7] T 2 +TiR 2 RiT 2 + TiR 2 RiR 2 RiT 2 + ... 

= 7\(1 +RiR 2 + (R { R 2 ) 2 + ...)T 2 = Tl l 2 u , (1.89) 

1 — K\I\2 
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or 



1 - T _ 1 - 7i 1 - T 2 
T ~ f t + T 2 ' 
which demonstrates the additivity of (l — T)/T. Thus, for N scatterers we obtain 



(1-90) 



1 - T 1 - 7i 

=AT -. 

r 7] 



(1.91) 



The resistance of the system is 

hi h h Ri 
R=^ = -+N~. (1.92) 

We again obtain the series resistance of N Landauer scatterers and contact resistance. 



1.3 Multi-channel scattering and transport 
1.3.1 S-matrix and scattering states 

Consider now the general multi-channel case, when the scattering is possible be- 
tween different modes. It is convenient to define separately left (L) and right (R), 
incoming (+) and outgoing (-) modes (Fig. 1.10). We assume that at z < Zl and z> Zr 
the leads have a constant cross-sections. To make the 5-matrix unitary, we introduce 
the normalization of incoming and outgoing modes, as was discussed in sec. 1.1.4. 
Thus, outside the scattering region we define 



^ (r) = vic fe(x,y)A " +e ^' Z<ZL 

VwbM = „L» 0/,,: v,V!-t„. c h 



\f2nhv Ln 

wLs( r ) = h L„ <j>R m (x,y)B m+e -< z 

y/27thV Rm 

&W = <$>R,n{x,y)B m -e ik '"< z 



\J2nhv, 



Z<ZL 

z>zr 
Z > Zr. 



(1-93) 
(1.94) 
(1-95) 
(1-96) 



'Rin 



A 



n+- 



A 



m- 



z 

Fig. 1.10. Multi-channel scattering. 
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Here n and m label the left and right transport channels respectively, (j>L n (x,y) 
and <pR m (x,y) are the left and right transverse eigenfunctions, Vlh =hk„/mi and 
VRm = hk' m /niR are the velocities, hk n — \j2mjjE —E„), hk' m — yj2mn(E — E m ), E n , 
E m being transverse eigenenergies. As in the single-channel case, we mark the left 
wave amplitudes with A± and right with B±. 
The initial incoming states become now 



V?(R)n E ( r ) = 



1 



^/2nhv L{R)n 



h(R)n(x,y)e 



+(-)ik n z 



(1.97) 



It is important to remember, that we change the normalization, but we must keep 
the density of states. Indeed, the matrix element (v^sI^e) acquires the additional 
multiplier 1 / (hv). Thus, we should change 



/§-/ 



dE 
hv 



dk J 
2% 



-hv = / dE 



This integration rule follows also from the normalization condition 



(1-98) 



(V l L(R)nE\V'L( R yE>) = J (VtJji^WJ Y% RWE ,(r)dxdydz = 8 nn ,8(E-E'). 

(1.99) 

Scattering matrix relates all incoming with all outgoing modes. In the single- 
channel case, considered previously, S is 2 x 2 matrix (1.31). In the general multi- 
channel case the 5-matrix is defined asN xN matrix 







(A 1+ \ 


A N L - 

Bi- 


= 8 


A N L + 

Bi+ 


\Bn r -J 




\ b n r + ) 



S21 S22 ■ 



■ Sin \ 

■ S2N 



( M + \ 

A N L A 

Bi+ 



(1.100) 



\5jvi Sn2 ■•■ Saw / \Bn r + ) 

where Nl and Nr are the numbers of left and right channels, N = Nl+Nr. 
We represent it in the block form, analogous to (1.31): 



S = 



(1.101) 



where the matrices f (Nl x Nl) and i (Nl x Nr) describe transmission and reflection 
of the states incoming from the left, the matrices ? (Nr x Nr) and V (Nr x Nl) de- 
scribe transmission and reflection of the states incoming from the right, Nl (Nr) is 
the number of the left (right) channels. 



The matrix t is called transmission matrix, T m 



are the probabilities of 



transmission from the left mode n into the right mode m, R n i n = \r n i t 
abilities of reflection from the left mode n into the left mode n', etc. 
The scattering marix is unitary 



are the prob- 



S^S = ss f = /, 



(1.102) 
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from which the conservation of total probability is clear 

(& § )nn = LIW| 2 + L|W| 2 = I- (1-103) 

n' m 

The other useful relations are 

r^r + Pt = r n r'+i n t' = /, (1.104) 

ff^ + f'f' t = r'r n + tP =1, ( 1 . 105) 

ft' + PP = t ,f r + r'H = 0, (1.106) 

rP + f r' f = t r 1 " + r'f /f =0. ( 1 . 107) 

From which follows 

Tr (Pt) = Tr (ff") = Tr (f' 1 ?') = Tr (f f f ) , (1.108) 

Tr (r + r) = Tr (rr f ) = Tr (f'V) = Tr (r'r' + ) . (1.109) 

In the time-reversal case the 5-matrix is symmetric 

S = S T , (1.110) 

thus the reflection matrix is also symmetric r nn i — r„i„ , and f ' = f r . 
In the external magnetic field more general conditions take place: 

r m ,{B) = r n , n {-B), (1.111) 

t nm {B)=t' mn {-B). (1.112) 

The simplest way to proceed is to use the so-called scattering states. Using in- 
coming and outgoing modes in the right and left leads (1.93)-(1.96), we can define 
the scattering states as 



/0 ^ <j> Rm (x,y)e ,k "- z , z > Z R . 
t, m V2nhv Rm 

(1.113) 

The physical sense of this state is quite transparent. It describes a particle moving 
from the left and splitting into reflected and transmitted parts. We established that 
it is important for Landauer transport, that only these states are populated from the 
left reservoir, so that one can accept the distribution of "left" scattering states to be 
equilibrium with the left electro-chemical potential. The "right" states, populated by 
the right reservoir, are defined as 



^ = ^ Jh ,(x,y)e- a ^ + £-^ = fa B ,(x,y)e tt ^, z > zr, 



^ /7 ^— hn(x,y)e- ik " z ,z<z L . 

(1.114) 

Now, using the left and right states(1.93)-(1.96) or the scattering states (1.113), 
(1.1 14), it is straightforward to obtain the Landauer formula. 
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1.3.2 Multi-channel Landauer formula 



Below in the section 1.3 we consider only time-reversal case, so that 

Rn'n = Vn'n\ 2 = \ r nn'\ 2 = Rim', (1.115) 

C- d-116) 



|2 _ u/ |2 _ 



Zero temperature conductance 

The heuristic calculation of the conductance is straightforward. The current from the 
left mode y^r) into the right mode yf Rrn (r) is determined by the transmission prob- 
ability T mn . The inverse current from \j/[ n (r) into \j/ Rm (r ) is determined by T' nm = T mn . 
Assuming that all left states y^r) are populated with the electrochemical potential 
of the left electrode, and right states W Rm {r) are populated with the electrochemical 
potential of the right electrode, we obtain the current of incoming left n and right m 
electrons as 

Imn = yT m n{flL-tlR), (1.117) 

h 

summing contributions from all incoming left modes and all incoming right modes 
we obtain the full current 

e N L N R 

1= t£ £ T mn (fL L -fL R ). (1.118) 

n n=\m=\ 

After the obvious mathematical transformation 

N L N R N L N R N L N R N L N R ^ 

Tmn = ^ ^ \tmn\ = 52 52 tmnt mn ~ 52 52 ) nm^ mn 
n=lm=\ n=\m=l n=lm=\ n=lm=\ 

N L 

= I(^L=Tr(f t f), (1.119) 

n=l 

we obtain the Landauer multichannel current and conductance 

2 

/ = ^Tr ft) (fe - ft,) = ^Tr (ft) V, ( 1 . 120) 

G=%-1r(iH). (1.121) 



General expression 

If we repeat the above summation procedure at finite temperature and voltage, we 
obtain a generalized form of the equation (1.63) 

1{ y) = zL r Tnm(E,V)\f L (E)-f R (E)}dE. (1.122) 

" nm 
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It is instructive, however, to consider an approach utilizing the scattering states. 
Let us calculate the current carrying by one scattering state ft„£(r) (1.113). Direct 
application of the current density expression 

ieTi 

j = ^ [V(r) VyA* (r) - y* (r) V y/(r)] (1. 123) 

gives the following answer for the z-component of the current density: 

f~ = sSs^^^'"-^-*^ (U24) 



And integrating over (x,y) and using the orthogonality property 

J j <$> Rm {x,y)<t> Rm ,(x,y)dxdy = 8 mm ,, (1.125) 

we get 



IimE = e£ \t mn \ 2 = I {?i) nn . (1.126) 

m " 

We can summarize over transport channels, and obtain the following general 
two-terminal Landauer formula for the current 

I(V)= e -J~J(E,V)[ME + eV)-ME)]dE, (1.127) 

with 

T(E,V) =E( ?+ 0„„ =Tr(P?) =Yjmn- (1-128) 

n nm 

T(E,V) is the effective transmission function for the particles with the energy E. The 
most important advantage of this formula is, that the transmission function can be 
calculated from the quantum scattering theory. Thus, the kinetic problem is reduced 
to the pure quantum mechanical problem of a single particle in a static potential. The 
formula (1.127) is the most general two-terminal formula. All other Landauer for- 
mulas can be obtained in the limiting cases from (1.127). For the finite-temperature 
conductance we have 



G= £ — ( /•:/•., ( (> ; 1 '' ) (1.129) 



o 2 {■<*• _ 

hLJ {E \ dE 
At zero temperature (—d fo/dE) = 8(E — Ep), thus 



G = —T(Ef). (1.130) 



In agreement with the previous results. 
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Transmission probabilities and transmission eigenvalues 

The physical sense of the transmission function (1.128) can be interpreted in two 
different ways. 

First, T n = (fit) is the transmission probability for the incoming state in the 
n-th channel. 

On the other hand, the NxN matrix Pi can be diagonalized. One can consider the 
eigenvalues ,% as the transmission coefficients in the basis of independent transport 
channels. 

The conductance has the same form in both cases because of the invariance of 
the trace 

G =tL T >i = jL^- (1.131) 



1.3.3 Derivation from the linear response theory 

In the previous sections we obtained the Landauer formula, valid at finite voltage, 
using some additional arguments about the populations of the states inside the junc- 
tion. Here we restrict ourselves by the linear response approach, which is valid only 
in the limit of small voltage, but is exact. 

We start from the Kubo formula for conductance (we follow Bruus and Flens- 
berg [8]): 

G= lim -Re H dte i ^ +i ^ t ( \l(t), 7(0)1 ) , (1.132) 

0)^0 CO Jo \ L J - / eq 

where [...]_ is the commutator, the matrix elements of the current operator / in the 
Hilbert space of single particle eigenfunctions (r) are 

which should be independent of z in the stationary case because of current conserva- 
tion, and can be evaluated at any cross-section. 

To proceed, we use the current operators in the second quantization representa- 
tion: 

t=Y,hvc\cv, (1-134) 

XV 

and choose the scattering states (1.113), (1.114) as eigenfunctions with X being 
equivalent to the set s,n,E, where s — L,R is the lead index, n is the number of 
channel in the lead s, and E is the energy. Straightforward calculation of the statistical 

average in equilibrium, using U\c X i \ = fo{E x )8 X \> and l xv = l* vv gives 



eq 



[/»,/(())] _) = £ vW (£i-£ *' 



■Y,\hx'\ 2 e i{E ^' ) \ME x )-fo(E^ 
XV 



eq 

(1.135) 
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Inserting this into (1.132) and integrating, we get 



G=lim h lmT 

m^o co fa co + ir] +Ex- Ey 



dfo 
dE 



h{E x )-f {E x ,) 
8{E X -E X '). 



(1.136) 



Now we use explicitly, that X =s,n,E, and integrating over the energy E' obtain 



dfo 
dE 



/co 
J^J^|/sB,s"l'l 
-OO / / 
ss nn 

= Kh ( £ [\I Ln .Ln'\ 2 + \I L n,Rn'\ 2 + l^Ln'^ + I 1 ! 

J — OO / 
WW' 



Rn.Rn' 



dfo 
dE 



dE. (1.137) 



The calculation of \I m y n >\ 2 for the scattering states (1.113), (1.114) is straightfor- 
ward, left or right parts of these expressions can be equivalently used, it is essentially 
the same, as we made to get (1.126). Thus, we have (here n,n' e L, m,m' € R) 

4,n' = ^ (f f f)„'„' = ft ( r ' tf 0n'n> 7 ™>« = - ^ (^Onn" l mrf = ~ ft (^) mm' ' 

(1.138) 

Using this expressions and the properties (1.105), it is easy to show that in the time- 
reversal case 

EEIW„<| 2 = (^) 2 Tr [{fff+^ff+ff^+^r] = 2 (0 2 Tr(ftf) , 

ss' nn 1 

(1.139) 

and finally 

dfo ' 



G = 



Tr(ftf) 



(1.140) 



We derived again the Landauer formula. 
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1.4 Multi-terminal systems 

1.4.1 Multi-terminal Landauer-Biittiker formula 

The scattering theory, developed in the previous section, can be applied also in the 
multi-terminal case (Fig. 1.11) exactly in the same way, as in the multi-channel case. 
The only difference from the considered before two-terminal system is that now 
there are several electrodes with independent electrochemical potentials and the same 
number of leads with independent transport channels. 

We use in this section the following notations: M is the number of terminals, 
N s (s = 1,...,M) is the number of transport channels in s lead, N = L. s ^.s is the full 
number of channels. 

The current through the terminal s, flowing in the direction to the scattering re- 
gion, is given by the following expression 



f°° dk f° 

/ V s (n,k z )f s {n,k z )— 1 + / V s {n,k z )Y\S m , pm \ f p (m,k z ) 
Jo 2n J-oo ^ 



2% 

(ii4i) 

here s, p are the terminal indices, n, m are the channel indices. The first term is the 
current of the electrons incoming from the s electrode, the second term is the current 
of electrons reflected and transmitted from other electrodes. Changing to the energy 
integration, we get finally 



J n.pm 



\ Z ]f P (E)dE. 



Let us introduce the transmission and reflection functions 



nm nn i 

by definition f ss = R ss . Then, the current gets the form 



I 2 

sn.sn' I j 



(N s -R„)ME)-Y i Ts P f p (E) 



dE, 



(1.142) 



(1.143) 



(1.144) 




Fig. 1.11. Multi-terminal system. 
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or, equivalently, 



N s f s (E)-Y,T sp f p (E) 



dE, 



or, equivalently 



I s = fl Lj_ oo i T Psfs(E) - T sp f p (E)\ dE. 
From the unitarity the following sum rule can be obtained 

y, = 52 t ps = n s . 

p p 

From time-reversal in the presence of a magnetic field H 

I^(H)=£f ps (-H). 



(1.145) 



(1.146) 



(1.147) 



(1.148) 



1.4.2 Biittiker conductance formalism 

Following Biittiker [9, 10], we formulate the elegant formulas, which allow to calcu- 
late the currents and voltages in the linear regime. We start from the current expres- 
sion (1.146). For small voltages 



i s =Z\- g p* v *- g *p v i 



p\ ■ 



(1.149) 
(1.150) 



In the time-reversal case these formulas take the usual multi-channel form. For 
example, the formulas (1.120) and (1.121) can be rewritten as 



2 2 



e 2 ^ e 2 



(1.151) 



G sp = - T sp = -Tr (t sp f sp ) . (1.152) 

In equilibrium I s = 0, but the potentials V s can have some constant value because 
of gauge invariance V s = Vq ^ 0. Thus, we obtain the sum rule 



p p 



(1.153) 



This shows, that the current can be equivalently written as 



/ s = Lg sp [v s -v p ]. 

p 



(1.154) 
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There is also the alternative representation 

p 

where the conductances G are defined by the following M xM matrix 



Ni-Ru -T l2 ■■■ -Tim 

— T21 N2-&22--- —T2M 

G' = 



(1.155) 



—Tmi 



■ ■ ■ Nm - Rmm 



(1.156) 



The conductances G and G' are related in an obvious way, and can be used equiva- 
lently. The sum rule for G is 



13 



ps 







g' 



LG'ps- 

p^s 



(1.157) 



Before discussing the applications of this current formula, we should introduce 
the notion of a voltage probe. The voltage probe is any terminal, in which the zero 
current is preserved, I s = 0. In spite of zero current, such terminal influences the 
current through other terminals. The physical reason is that the electrons can enter 
the electrode and return back, but after phase randomization. Thus, the voltage probe 
can be used not only for voltage measurement, but also for introducing a decoherence 
into the system. 



3-terminal conductance 

Consider as an example, a three-terminal junction (Fig. 1.12). Let us consider the ter- 
minal 3 as a voltage probe: I3 = 0, h = —I2 = I, Vi — V2 = V. We have the following 
equations 



h = G31 (V3 - Vi) + G 32 (V3 - V 2 ) = 0, 
GnV + G13 (Vi - V 3 ) = G21V - G 23 (V 2 -V 3 ). 



We can also choose V\ ~V,V2~ 0, then 



V 3 



G31 



G31 +G32 



V. 



(1.158) 
(1.159) 



(1.160) 



This voltage is not small in the limit G3i,G32 <C G12, when the third electrode is 
used for unperturbed voltage measurement. The current is 



/ = 



G12 + 



G13G32 



G31 +G32 



V. 



(1.161) 



it is not perturbed when G3i,G32 — > 0, but at finite coupling to the third electrode, 
the current between electrodes 1 and 2 is changed. 
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4-terminal conductance 

Consider now the 4-terminal system (Fig. 1.12) and assume that electrodes 3 and 4 
are weakly coupled voltage probes. Ih this case the current / = I\ = —h ~ GnV . 
What is the voltage V' — V3 — V4? Neglecting G34 and G43 as the next-order small 
terms, we obtain 

yl= (highly (1162) 

(j31 +tJ32 

and 



/ _ e 2 R 

We obtained the first Landauer formula. 



(1.163) 



2 1 



Fig. 1.12. 3 and 4 terminal junctions. 



1 .4 Multi-terminal systems 

1.4.3 Edge states in quantum Hall effect 
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Problems 

1.1 Inter-mode scattering in an adiabatic waveguide 

Find the nonadiabatic corrections A n to the equation (1.7) 

ft2 32 



{ ~ Jz 2 + E " ^ } WnE ^ = E WnE ^ + A " ' 



1.2 Transmission though a symmetric rectangular barrier 

Consider the rectangular barrier 

{0, z < —a, 
U , -a<z<a, 
0, z>a. 

The wave function is given by 

( A + e ikz +A-e- ikz ,z<-a 
\lf(z) = < Cie- Kz +C 2 e K \ -a<z<a 
\B-e ikz +B+e- ikz ,z>a 



where 



y/2mE z 



Show that 



T(E Z < U ) 



T(E Z > U ) = 



ft 1 

^/2m(U -E z ) 

ft 

4k 2 K 2 



(k 2 + K 2 ) 2 smh 2 (aK)+4k 2 K 2, 
Ak 2 K 2 

{k 2 -K 2 ) 2 sva 2 {aK)+Ak 2 K 2 ' 



1.3 Transmission though a 5-barrier with asymmetric leads 

1.4 Transmission though an asymmetric rectangular barrier 
Other problems 



